Quantum Trajectory method for the Quantum Zeno and anti-Zeno effects 
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We perform stochastic simulations of the quantum Zeno and anti-Zeno effects for two level system 
and for the decaying one. Instead of simple projection postulate approach, a more realistic model 
of a detector interacting with the environment is used. The influence of the environment is taken 
into account using the quantum trajectory method. The simulation of the measurement for a 
single system exhibits the probabilistic behavior showing the collapse of the wave-packet. When a 
large ensemble is analysed using the quantum trajectory method, the results are the same as those 
produced using the density matrix method. The results of numerical calculations are compared with 
the analytical expressions for the decay rate of the measured system and a good agreement is found. 
Since the analytical expressions depend on the duration of the measurement only, the agreement 
with the numerical calculations shows that other parameters of the model are not important. 

PACS numbers: 03.65.Xp, 42.50.Lc, 03.65.Ta 



IT) . I. INTRODUCTION 

The quantum Zeno effect has attracted much attention. The effect is caused by the influence of the measurement on 
| the evolution of a quantum system. This effect is related to the non-exponential survival probability of the quantum 
£N) ■ system. 

The exponential survival law is known to be an excellent phenomenological fit to unstable phenomena. However, 
from quantum mechanics it follows that the survival probability is not exponential for short and long times. Short- 
time behavior of the survival probability is not exponential but quadratic 0. The deviation from the exponential 
decay has been confirmed by Wilkinson et al. 0. This result, when combined with the frequent measurements, 
leads to what is known as the quantum Zeno effect Q . Nowadays there are a number of experiments which claim to 
have verified the quantum Zeno effect and some others are planned 0, 0, IS • It was also predicted that frequent 
measurements could accelerate the decay process |cl IgL Hoi fill IT2I fli! . This is the so-called quantum anti-Zeno effect. 
Both effects were first observed in an atomic tunneling experiment |l4j. 

The states of the system need not to be frozen: in the general situation the coherent evolution of the system can take 
place in dynamically generated quantum Zeno subspaces [l5j . The projective measurements used in the description 
q-i of the quantum Zeno effect can be replaced by another quantum system interacting strongly with the principal one 

. ^ \ Interaction with the measuring device is one of many possible interactions of the system with an external environ- 
ment. It is known that not only measurements cause consequences similar to the quantum Zeno effect on the system's 
evolution The experiment of Itano et al. Q has been explained in Refs. [2^, H3, using the dynamical 

description, without using the concept of the measurements. It was shown that the quantum Zeno effect follows from 
the quantum theory of irreversible processes, as well. Therefore, the quantum Zeno and anti-Zeno effects can be 
considered as more general phenomena. However, in this paper we will consider only one particular phenomenon, i.e., 
the effects caused by the measurements. 

Quantum Zeno effect can have practical significance in quantum computing. The use of the quantum Zeno effect for 
correcting the errors in quantum computers was first suggested by Zurek |2jj. A number of quantum codes utilising 
the error prevention that occurs in the Zeno limit have been proposed |2a . 1271 1281 129L l30| . 

In the analysis of the quantum Zeno effect the projection postulate is not sufficient. The measurement should be 
described more fully, including the detector into the description. In the description of the quantum measurement 
process, the environmentally induced decoherence plays a very important role [3lll32l l33l l34l I35ll3rll3^ | . Therefore, 
in order to correctly describe the measurement process one should include into the description the interaction of 
the detector with the environment. In this paper we describe the evolution of the detector interacting with the 
environment using the quantum jump model developed by Carmichael [38|. 
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The density matrix analysis assumes that the experiment is performed on a l arge number of systems. An alternative 
to the density matrix analysis are stochastic simulation methods [38l l39l I40L fill . Various stochastic simulation 
methods describe quantum trajectories for the states of the system subjected to random quantum jumps. Using 
stochastic methods one can examine the behaviour of individual trajectories, therefore such methods provide the 
description of the experiment on a single system in a more direct way. The results for the ensemble are obtained by 
repeating the stochastic simulations several times and calculating the average. 

Stochastic simulations of the quantum Zeno effect experiment were performed in Ref. [4^ | . In present paper we use 
the quantum jump method to describe the evolution of frequently measured systems and to compare the numerical 
results with the analytically obtained decay rates. 

The paper is organized as follows. In Sec. [n] we present the description of the measurement. The model of the 
detector is presented in Sec. IIIII In Sec. |W] we present the quantum jump method briefly. In Sec. E] evolution of 
the detector is calculated using the quantum jump method. Evolution of frequently the measured two level system 
is investigated in Sec. IVII In Sec. IVHI we present a numerical model of the decaying system. Using this model in 
Sec. lVIlil we calculate the evolution of the frequently measured decaying system. Section llXl summarizes our findings. 



II. DESCRIPTION OF THE MEASURED SYSTEM 



We consider a system that consists of two parts: A and F. The system A is interacting with the detector, i.e., it 
is measured. We assume that the system A has the discrete energy spectrum. The Hamiltonian of this part is Ha- 
The other part of the system is represented by Hamiltonian Hf. Hamiltonian Hf commutes with Ha- The operator 
V causes the jumps between different energy levels of Ha- Therefore, the full Hamiltonian of the system is of the 
form Hs = Ha + Hp + V. The example of such a system is an atom with the Hamiltonian Ha interacting with the 
electromagnetic field, represented by Hp, while the interaction between the atom and the field is V. 

In this article we consider the system A that has two levels: ground \g) and excited |e). We will measure whether 
the system is in the ground state. The measurement is performed by coupling the system with the detector. The full 
Hamiltonian of the system and the detector equals to 

H = Hs+H D +H T , (1) 

where Hp, is the Hamiltonian of the detector and Hj represents the interaction between the detector and the measured 
system, described by the Hamiltonian Ha- We can choose the basis \na) = \n) \a) common for the operators Ha 
and Hp, 

H A \n) = huj n \n), (2) 
Hp\a) = Hu> a \a). (3) 

Here tvjj n and hu a are energies of the systems A and F, respectively. 

The initial density matrix of the system is /3s (0). The initial density matrix of the detector is pp(0). Before 
the measurement the measured system and the detector are uncorrelated, therefore, the full density matrix of the 
measured system and the detector is p(0) — ps(Q) ® Ad(0). 

When the interaction of the detector with the environment is taken into account, the evolution of the measured 
system and the detector cannot be described by a unitary operator. More general description of the evolution, allowing 
to include the interaction with the environment, can be given using the superoperators. Therefore, we will assume 
that the evolution of the measured system and the detector is given by the superoperator S(t). 



A. Measurement of the unperturbed system 

In this section we investigate the measurement of the unperturbed system, i.e., the case when V = 0. 

We assume that the Markovian approximation is valid i.e., the evolution of the measured system and the detector 
depends only on their state at the present time. The master equation for the full density matrix of the detector and 
the measured system is 

£p(t) = hH A) p(t)} + \[Hj,p(t)} + \[H D ,p(t)]+£ D p(t), (4) 
at in in in 

where the superoperator Hp, accounts for the interaction of the detector with the environment. We assume that the 
measurement of the unperturbed system is a quantum non-demolition measurement |44l l45l EH |47| . The measurement 
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of the unperturbed system does not change the state of the measured system when initially the system is in an 
eigenstate of the Hamiltonian H^. This can be if [Ha, Hj] = 0. 

We introduce the superoperator £„ )m acting only on the density matrix of the detector 

Cn, m p = ^r((n\Hi\n)p - p(m\Hi\m)) + ^r[H D ,p] + C D p (5) 
in in 

and the superoperator <S„ im (t) obeying the equation 

d 

with the initial condition <S n , m (0) = 1. Then the full density matrix of the detector and the measured system after 
the measurement is 

P (tm) = S(t m ) P (0) = \n)(ps)n, m (0y^™(m\ <g> 5„, m (rM)^(0), (7) 
where tm is the duration of the measurement and 

Wm : n = U) m — Ul n (8) 

with uj n defined by Eq. J3J). From Eq. JJJ it follows that the non-diagonal matrix elements of the density matrix of 
the system after the measurement (ps)n,m( T M) = {ps)n,m{fy e%ulm nTM are multiplied by the quantity 

F, hm (T M ) = Ti'{5„, m (r M )/5z)(0)}. (9) 

Since after the measurement the non-diagonal matrix elements of the density matrix of the measured system should 
become small (they must vanish in the case of an ideal measurement), F ni . m (TM) must be also small when n ^ ni. 



III. THE DETECTOR 



We take an atom with two energy levels, the excited level \a) and the ground level \b), as the detector. The 
Hamiltonian of the detecting atom is 

H D = ^a z . (10) 

Here MId defines the separation between levels \a) and \b), a x , a y , a z are Pauli matrices. The interaction Hamiltonian 
Hi we take as 

H I = h\\g){g\(a + + a-), (11) 

where a± = \(o x ± i&y)- The parameter A describes the strength of the coupling with the detector. The detecting 
atom interacts with the electromagnetic field. The interaction of the atom with the field is described by the term 

r 

CdPd = --(a+cf-pD - 2a-p D a + + p D a + a-), (12) 

where T is the atomic decay rate. 

At the equilibrium, when there is no interaction with the measured system, the density matrix of the detector is 
p D (0) = \b){b\. 



A. Duration of measurement 



We can estimate the characteristic duration of one measurement tm as the time during which the non-diagonal 
matrix elements of the measured system become negligible. Therefore, in order to estimate the duration of the 
measurement tm, we need to calculate the quantity 



F e:g (t)=Tv{S e Jt) PD (0)}. 
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We will solve the equation 



d 1 

-qJd = T^[H D ,PD]+£e,gPD. (13) 

For the matrix elements of the density matrix of the detector we have the equations 

d 

TTTPaa = l^Pab ~ ^ Pan, (14) 

at 
d 

-Q t Pbb = l^Pba+^Paa, (15) 
d 1 

Q^Pab = -i^DPab + l^Paa - ^ Pab, (16) 
8 1 

— p ba = iVl D p ba + i\p bb - -Tp ba . (17) 

with the initial conditions p a b(0) = Pba(0) — Paa(0) — 0, pbb(0) = 1- 

Atom can act as an effective detector when the decay rate T of the excited state \a) is large. In such a situation we 
can obtain approximate solution assuming that p ba and p ab are small and p bb changes slowly. Then the approximate 
equation for the matrix element p ba is 

d 1 

foPba = «Ap b& (0) - -Tp ba (18) 

with the solution 

p 6a (t) = 2^(l-e-^). (19) 
Substituting this solution into the equation for the matrix element p bb we get 

J^6 = -2^(1 -e-i"). (20) 
Taking the term linear in t we obtain the solution 

Pbb (t) » 1 - 2yt » e" 2 ^* . (21) 

Since p aa {t) = 0, using Eq. we have 

F e , g (t) = pbb(t) ~ exp 

where 



t 

TM 



tm = ^ (22) 

is the characteristic duration of the measurement. This estimate is justified comparing with the exact solution of the 
equations. We get that the measurement duration is shorter for bigger coupling strength A. 



IV. STOCHASTIC METHODS 



The density matrix approach describes the evolution of a large ensemble of independent systems. The observed 
signal allows us to generate an inferred quantum evolution conditioned by a particular observed record |38j . This 
gives basis of the quantum jump models. In such models the quantum trajectory is calculated by integrating the 
time-dependent Schrodinger equation using a non-Hermitian effective Hamiltonian. Incoherent processes such as 
spontaneous emission are incorporated as random quantum jumps that cause a collapse of the wave function to a 
single state. Averaging over many realizations of the trajectory reproduces the ensemble results. 
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The theory of quantum trajectories has been developed by many authors IS m E3, H El US i3. Quantum 
trajectories were used to model continuously monitored open systems [3g, I3 |4£| , in the numerical calculations for 
the study of dissipative processes [H 03] , and in relation to quantum measurement theory 0, ■ 

We assume that the Markovian approximation is valid. The dynamics of the total system consisting of the measured 
system and the detector is described by a master equation 

= (23) 
where M. is the superoperator describing the time evolution. The superoperator M. can be separated into two parts 



M = £ + J. (24) 

The part J is interpreted as describing quantum jumps, C describes the jump-free evolution. After a short time 
interval At the density matrix is 

p(t + At) = p(t) + Cp(t)At + Jp(t)At. (25) 
Since Eq. (|23[l should preserve the trace of the density matrix, we have the equality 

Tr{£p(t)} + Tr{Jp(t)} = 0. (26) 
Using Eq. (126[) equation l|25|) can be rewritten in the form 

* + - f&X* <' - ^*»"> + 4jm m]M - <27) 

This equation can be interpreted in the following way: during the time interval At two possibilities can occur. Either 
after time At the density matrix is equal to conditional density matrix 



with the probability 



or to the density matrix 



Pjump(t) = Tr{Jp{t)}At (29) 
pit) + CAtpjt) 

P„o- Jump (t + At) = rp —— (30) 

with the probability 1 — pj UII1 p(t). Thus the equation ^jfy can be replaced by the stochastic process. 
Here we assume that the superoperators C and J have the form 

Cp = ±(H eS p-pHl s ), (31) 
Jp = Cp&. (32) 

The operators H e g and C are non-Hermitian in general. If the superoperators C and J have the form given in 
Eqs. <|3 ljl. (|32[) and the density matrix at the time t factorizes as p(t) = \^(t)) (&{t)\ then after time interval At the 
density matrices p iump (t + At) and j5 no _j ump (t + At) factorize also. Therefore, equation for density matrix (|27ll can 
be replaced by the corresponding equation for the state vectors. The state vector after time At in which a jump is 
recorded is given by 

|* jump (t + At)) = 1 =C\*(t)). (33) 

^<*(t)|£t(?|*(t)) 

The probability of a jump occurring in the time interval At is 

Pjump(t) = (*(t)|C t (7|*(t))At. (34) 
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If no jump occurs, the state vector evolves according to the non-Hermitian Hamiltonian _ff c ff, 

|¥ jump (t + At)) = 1 (l - i-HeffAt ) |*(/ )). ( V>) 



Numerical simulation takes place over discrete time with time step At. When the wavefunction |\&(£ n )) is given, 
the wavefunction |*(i„ + i)) is determined by the following algorithm |38| : 

1. evaluate the collapse probability p jnmv (t n ) according to Eq. (|34[1 

2. generate a random number r n distributed uniformly on the interval [0, 1] 

3. compare pj ump (i n ) with r n and calculate \^ c (t n+ i)) according to the rule 

|*(*n+0) ~ exp(-~H BS At \ |¥(t»)), Pjump(tn) > T„. 

We can approximate the second case as 

|*(*n+l)) ~ (l - At ) !*(*«)>• 

V. STOCHASTIC SIMULATION OF THE DETECTOR 

At first we consider the measurement of the unperturbed system and take the perturbation V = 0. The measured 
system is an atom with the states \g) and |e). The Hamiltonian of the measured atom is 

H A = tkJ A \e)(e\, (36) 

where Hlua is the energy of the excited level. 

The stochastic methods described in Sec. IIVI were used to perform the numerical simulations of the measurement 
process. Using the equation l|12[) we take the operator C in Eq. I|32|) describing jumps in the form 

C = Vf er_ (37) 

and the effective Hamiltonian in Eq. (|31|) as 

r 

iJoff = Ha + Hjj + Hj — ih-a+a- . (38) 

The wavefunction of the measured system and the detector is expressed in the basis of the eigenfunctions of the 
Hamiltonians of the measured system and the detector 

|*) = c ea \e)\a) + c eb \e)\b) + c ga \g)\a) +c gb \g)\b). (39) 

The effective Hamiltonian produces the following equations for the coefficients of the wave function I*) 

C ea = -i + - i^j C ea , (40) 

Ceb = ~i (uA - -^-J Ceb, (41) 



o_d r 

2 2 - 



-i\c gb - i— c ga - -c ga , (42) 



-i\c ga + i^Y-c gb . (43) 
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FIG. 1: Typical quantum trajectories of the detector. Figure shows the probability p aa of being in the excited level of the 
detector. The solid line corresponds to the case when the measured system collapses to the ground state and the dashed line 
corresponds to the case when the measured system collapses to the excited state. The parameters used for numerical calculation 
are At = 0.1, F = 10, A = 1, and tl D = I. 

Equations (|40|) - (|43|l are used in the numerical simulations to describe the evolution between the jumps. After the 
jump in the detecting atom the unnormalized wavefunction is 

C\y)=VT(c ea \e)\b)+c ga \g)\b)). (44) 

The jump occurs with the probability Pj um p obtained from Eq. I|34|) . 

D . - YAt |Cea| 2 + \Cgg\ 2 , . 

For numerical simulation we take the measured system in an initial superposition state + 1 9))- The typical 

quantum trajectories of the detector are shown in Fig. ^ There are two kinds of trajectories corresponding to the 
collapse of the measured system to the excited or the ground states. The trajectories corresponding to the collapse 
of the measured system to the ground state show the repeated jumps. The mean interval between jumps, obtained 
from the numerical simulation is of the same order of magnitude as tm = 5 according to Eq. I|22|l . After averaging 
over the realizations the probability for the detector to be in the excited state is shown in Fig. [21 The figure shows 
that this probability reaches the stationary value. The time dependency of the non-diagonal matrix elements of the 
density matrix of the measured system is shown in Fig. [3] The figure shows a good agreement between the results of 
numerical calculations and the exponential decay with the characteristic time estimated from Eq. 11' I'D . 



VI. FREQUENTLY MEASURED PERTURBED TWO LEVEL SYSTEM 



We consider an atom interacting with the classical external electromagnetic field as the measured system. Interaction 
of the atom with the field is described by the operator 

V = -m R (\e)(g\ + \g)(e\)cosnt, (46) 

where f2 is the frequency of the field and is the Rabi frequency. In the interaction representation and using the 
rotating-wave approximation the perturbation V is 

V(t) = -n^(e iA "*|e)( 5 | + e -^| 5 )( e |), (47) 
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FIG. 2: Probability for the detector to be in the excited state, after performing an ensemble average over 1000 trajectories. 
The parameters used are the same as in Fig. 



where 

Aw = w A -0 (48) 

is the detuning. Here fku^ = hw e — huj g is the energy difference between the excited and the ground levels of the 
measured atom. 

If the measurements are not performed, the atom exhibits Rabi oscillations with the frequency Qr. If the measured 
atom is initially in the state \g), the time dependence of the coefficient c g of the wavefunction |\&) = c e |e) + c g \g) is 



,(t) = e -^ itA " ^cos Q*y&To| 



Aw /l 



In particular, if the detuning Auj is zero, we have 



c g (t)=c g (0) cos (^*). (50) 

When the measured atom interacts with the detector, we take the wavefunction of the measured system and of the 
detector as in Eq. I|39|) . In the interaction representation the equations for the coefficients, when the evolution is 
governed by the effective Hamiltonian H c h, defined by Eq. (|31(l . are 

(51) 
(52) 



Cea 


■^■R it Auj n 

— ^ n C Cqa 

2 


.fir. 


r 

- ~c ea 


Ceb 


it Auj . 

= * — e c sf> + 


.fiu 

i C e b 

2 




c ga 


o 


- iAc g b - 


,fifl 

* 2 Csa 


Cg b 


•^fl -itAw 

= t— e c e6 


— iXc ga 4 


n D 

l — Cgb 



Cga (53) 

(54) 



The evolution of the measured atom significantly differs from the Rabi oscillations. We are interested in the case 
when the duration of the measurement tm is much shorter than the period of Rabi oscillations 27r/£7^. In such a 
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300 



FIG. 3: Non-diagonal matrix elements of the density matrix of the measured system. The solid line corresponds to the 
numerical calculations and the dashed line corresponds to the exponential decay with the characteristic time given by Eq. t!Hl . 
The parameters used are the same as in Fig. Q 



situation the non diagonal matrix elements of the density matrix of the measured system remain small and the time 
evolution of the diagonal matrix elements can be approximately described by the rate equations 



fjj-Psa ~ F e ~*gPee(t) T g^ePgg{t), 

If the measured atom is initially in the state \g), the solution of Eqs. I|55|) and (|56|l is 



PggM = r , r = 9 i 1 + e 

1 e — >g 1 ^ g — >e ^ 

We can estimate the rates r e _> g and r g ^ e using equations from Ref. [5^ |. i.e 



(55) 
(56) 

(57) 



r 



e^g 



2tt 
2tt 



G(u)P eg (u))duj, 
G{u>)P ge (uj)dui, 



(58) 
(59) 



where 



Pge{u) 



1 f°° 

-Re / F eg (T)e j(w -^ )T dr 







-Re 

7T 



F ge (T)e t(ul+UJA)T dT 



(60) 
(61) 



and 



G(uj) = ( ) [<5(w - + Aw) + <5(w + u) A - Auj)) . 



(62) 
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FIG. 4: Typical quantum trajectory of the measured two level system (solid line). Figure shows the probability p gg for the 
measured atom to be in the ground level. The dashed line shows the Rabi oscillations in the measurement-free evolution. The 
detuning Auj is zero. The parameters used for the numerical calculation are At = 0.1, V = 10, A = 1, Qd = 1, and Qr = 0.1. 



Here the expression for G{ui) is derived using Eq. I|47|) . In contrast to Ref. [53] in the expression for P(lu) we extended 
the range of the integration to the infinity since F eg (T) naturally limits the duration of the measurement. Expressions, 
analogous to l|58|l . were obtained in Refs. as well. 

Using Eqs. (|58|) - (|62[1 we can estimate the transition rates as 

r„ s - r„ - ^ Ref ^ Me — * . f jf - f TTT ^ F (63) 

Here we used the expression exp(— t/tm) for Fegi^)- When the detuning Aw is zero the transition rates are 

r^,«r^«^H. (64) 

The transition rates are smaller for the shorter measurements. This is a manifestation of the quantum Zeno effect. 

For numerical simulation we take the measured system in an initial state \g). The typical quantum trajectory of the 
measured system is shown in Fig.0J The behaviour of the measured system strongly differs from the measurement- free 
evolution. The measurement-free system oscillates with the Rabi frequency, while the measured system stays in one 
of the levels and suddenly jumps to the other. The probability for the measured atom to be in the ground state 
calculated after averaging over the realizations is shown in Fig. The figure shows that this probability exhibits 
almost the exponential decay and after some time reaches the stationary value close to 1/2. The figure shows a good 
agreement between the results of the numerical calculations and the estimate (|57f) . 

When the detuning Aw is not zero, the frequently measured two level system can exhibit the anti-Zeno effect. This 
is pointed out in Ref. • For the case of nonzero detuning the probability that the atom is in the ground state is 
shown in Fig. [fj] The figure shows that the probability for the atom to be in the initial (ground) state is smaller, 
and, consequently, to be in the excited state is greater when the atom is measured. This is the manifestation of the 
quantum anti-Zeno effect in the two level system. 
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200 400 600 800 1000 

t/At 

FIG. 5: Probability for the measured atom to be in the ground level, after performing an ensemble average over 1000 trajectories. 
The solid line shows the results of the numerical calculations, the dashed line shows the approximation according to Eq. 1571 . 
The parameters used are the same as in Fig. |1] 




t/At 

FIG. 6: Probability for the atom with the nonzero detuning to be in the ground level. The solid line shows the results of the 
numerical simulation. The dotted line shows the approximation according to Eq. 15711 . using the decay rate from Eq. 1631 . The 
dashed line shows the evolution of a not measured system. The parameters used for numerical calculation are At — 0.001, 
r = 10, A = 1, Q.d = 1, Mr = 0.1, and Alu = 0.2. 
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VII. DECAYING SYSTEM 

We model the decaying system as two level system A interacting with the reservoir F consisting of many levels. 
The full Hamiltonian of the system is 

H = H A + H F + V, (65) 

where 

H A = hw A \e)(e\ (66) 

is the Hamiltonian of the two level system, 

H F =J2tuok\k)(k\. (67) 

k 

is the Hamiltonian of the reservoir, and 

V = h^2(g(k)\e)(k\+g(kr\k)(e\) (68) 
fe 

describes the interaction of the system with the reservoir, with g{k) being the strength of the interaction with reservoir 
mode k. In the interaction representation the perturbation V has the form 

V(t) = e i{H A +H F )ty e -i{H A +H F )t = %Y j {g{k)e i{ - UJA - u ^ t \e){k\ + g(k)* e - l(u,A - u)k)t \k){e\). 

k 

The wavefunction of the system A + F in the interaction representation is expressed as 

|*)= Ce (*)|e)|0)+X;pfcb)|A)- (69) 

k 

One can then obtain from the Schrodinger equation the following equations for the coefficients 

c e = -^jWe'^^^ct, (70) 

k 

t k = -igikfe-^-^Cz. (71) 
The initial condition is \^) — |e)|0). Formally integrating the equation (|71fl wc obtain the expression 

c k = -ig(k)* f e -*(«-A-«-*)f ' Ce (t')dt' '. (72) 



Inserting Eq. I|72|l into Eq. I|7l)|l , we obtain the exact integro-differential equation 

j f c e = -^dt'J2\9(k)\ 2 e i ^-^ t - t '^c e (t'). (73) 
The sum over k may be replaced by an integral 

— > / duj k p{uJk) 

k J 

with p(u>k) being the density of states in the reservoir. The integration in Eq. I|73(l can be carried out in the Weisskopf- 
Wigner approximation. We get the equation 

d r (0) 

5 Ce = _i£l^ Ce , (74) 
dt 2 v ; 

where the decay rate Fe-*g is given by the Fermi's Golden Rule: 

r(°i ff = 2^(^)| ff ( WA )| 2 . (75) 
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In order to observe the quantum anti-Zeno effect one needs to have sufficiently big derivative of the quantity 
p(uj)\g(uj)\ 2 . In such a case the decay rate given by Fermi Golden Rule l|75|l is no longer valid. The corrected decay 
rate may be obtained solving Eq. 1)73(1 by the Laplace transform method |54j . The the Laplace transform of the 
solution of Eq. (|73|1 is 



where H(z) is the resolvent function 



U*) = W y (76) 

H(z) = z+[ ^' g( ^V (77) 
J z + i{uu-uu A ) 

In the numerical calculations we take the frequencies of the reservoir ui distributed in the region [oja — A, oja + A] with 
the constant spacing Au>. Therefore, the density of states is constant p(uj) — 1/ Aui = po- The simplest choice of the 
interaction strength g(uj) is to make it linearly dependent on to, 

g(u) = g (l + j(u - wa)) , (78) 

where a is dimensionless parameter. Using Eq. I(78|l one obtains the expression for the resolvent 

H{z) = z + irp g% 1 - 2 arctan (|) + (a 2 | - i2a) - | + ~ arctan (|)) 

The real part of z at which the resolvent TL{z) is equal to zero gives the decay rate. Expanding the resolvent into the 
series of powers of A" 1 and keeping only the first-order terms we obtain the decay rate 

r« 9 -rn^i-^( 5a 2 -i)). (79) 

We solve Eqs. (|7U|I and (|71fl numerically, replacing them with discretized versions with the time step At. For 
calculations we used TV = 1000 levels in the reservoir. The numerical results for constant interaction strength are 
g(k) = go presented in Fig. The figure shows a good agreement between the numerical results and the exponential 
law according to the Fermi's Golden Rule at intermediate times. At very short time the occupation of the excited 
level exhibits quadratic behaviour, which, for the repeated frequent measurements, may result in the quantum Zeno 
effect. 

Numerical results in the case when the interaction with the reservoir modes is described by Eq. (|78|l with nonzero 
parameter a are presented in Fig. [S] The figure shows good agreement between the numerical results and the 
exponential decay with the decay rate given by Eq. I|79|l at intermediate times. For very short time we observe 
the acceleration of the decay due to the interaction with the reservoir. This acceleration for the repeated frequent 
measurements results in the quantum anti-Zeno effect. 



VIII. MEASUREMENT OF THE DECAYING SYSTEM 



In this section we consider the decaying system, described in Sec. IVIII and interacting with the detector. The 
wavefunction of the measured system and the detector, when the detector interacts with the decaying system in the 
ground state only, we take in the form 

|* c > = c ea \e)\0)\a) +cei|e)|0)|6) + ^(c fea | 5 )|fc)|a) + c kb \g)\k)\b)). (80) 

k 

The equations for the coefficients, when the evolution is governed by the effective Hamiltonian if e ff, are 



-i 9{k)e l ^ A -^ t c ka - t^c ea - | c eQ , 

k 



c eb = -iJ29(k)e i(uA - Uk)t c kb + i^c eb) (82) 



k 



Cka 
Ckb 



Mk)*e^ A ^ k)t c ea - i\c kb - i^-c ka - -c ka , (83) 
-i 5 (fc)*e-^-^* Ceb - i\c ka + i^c kb . (84) 



14 




FIG. 7: Time dependence of the occupation of the exited level of the decaying system. Solid line shows the results of the 
numerical calculation, dashed line shows the exponential decay according to the Fermi's Golden Rule. The parameters used 
for the numerical calculation are At = 0.1, Aui = 0.001, A = 0.5, and go = 0.001262. For the parameters used the decay rate 
is rglg = 0.01. 




FIG. 8: Time dependence of the occupation of the exited level of the decaying system when the interaction with the reservoir 
modes is described by Eq. 1781 . Solid line shows results of numerical calculation, dashed line shows exponential decay using 
the decay rate from Eq. (1791 . dotted line shows exponential decay according to Fermi's Golden Rule 17511 . In the calculations 
we used a = 2, other parameters are the same as in Fig. Q 
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After the jump in the detecting atom the unnormalized wavefunction becomes 

C|* c ) = \/f(c ea |e)|0)|6) +^2c ka \g)\k)\b)). (85) 

k 

According to Ref . 52) , the decay rate of the measured system is given by expression l|58|l with 

G( W ) = p( W )| ff ( W )| 2 (86) 

and 

P(w) = -Re/ F eg (T)e l ^- UA ^dT. (87) 

7T Jo 

Using F eg (r) = exp(— t/tm) we obtain 



P(w) 



1 TM 



7T 1 + (W - UJ A ) 2 Tl 



' M 

In order to obtain the quantum Zeno effect we take G(u>) as a constant 

fc2 2 

G(uj) = —™-, u> A -A<u><w A + A. (89) 
Alo 

Here Aui is the spacing between the modes of the reservoir. Using Eq. (|58|l we get the decay rate of the measured 
decaying system 

r e ^ 9 = TfX q - arctan(AT M ). (90) 



When Atm is big, we obtain the expression 



2 1 



r - = I< -( 1 -^ + -J (91) 

by expanding Eq. (|90|l into series of the powers of (Atm-) -1 . The second term in Eq. (|91|l shows that the decay rate 
decreases with the decreasing duration of the measurement tm- This is the manifestation of the quantum Zeno effect. 

The results of the numerical simulation are presented in Figs. El and 1101 Typical quantum trajectory of the 
measured decaying system is shown in Fig. EI This trajectory is compatible with the intuitive quantum jump picture: 
the system stays in the excited state for some time and then suddenly jumps to the ground state. The probability that 
the measured system stays in the excited state is presented in Fig. 1101 Figure shows a good agreement between the 
numerical simulation and the exponential law approximation with the exponent given in Eq. (|90|) . Also the quantum 
Zeno effect is apparent. 

When the detector interacts with the excited state of the decaying system the interaction term is 

H x = h\\e){e\(a+ +ct_) (92) 

and the quantum trajectories are different. Typical quantum trajectory is shown in Fig. 1111 This difference can be 
explained in the following way: when the detector interacts with the ground state, the interaction effectively begins 
only after some time, when the probability to find the system in the ground state is sufficiently big. This explains the 
absence of the collapses at short times in Fig. EJ Then, the measurement result after the collapse most likely will be 
that the system is found in the ground state. When the detector interacts with the excited state of the system, the 
interaction starts immediately and soon after that the most probable result of the measurement is that the measured 
system is in the excited state. It should be noted, that the averaged evolution shown in Fig. ^| does not depend on 
the state the detector is interacting. 

The model used for the decaying system when g{ui) — const does not allow to obtain the quantum anti-Zeno effect, 
since the conditions for the quantum anti-Zeno effect, presented in Ref. [9j, are not satisfied. In order to obtain the 
quantum ant-Zeno effect we use the interaction with the reservoir modes described by Eq. J7HJ. In such a case we 
have 

j-2 2 2 

G(w) = (l + j(w-w A yj , lo a -A<lu<lu a + A. (93) 
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FIG. 9: Typical quantum trajectory (solid line) of the measured decaying system when the detector interacts with the ground 
state of the measured system. Figure shows the probability p ee that the measured atom is in the excited level. The dashed line 
shows the exponential decay according to the Fermi's Golden Rule in the measurement-free evolution. The parameters used for 
the numerical calculation are At = 0.1, V = 10, A = 1, £Id = 1, Alj = 0.001, A = 0.5, and go = 0.001262. For the parameters 
used the decay rate is T e ->g = 0.01. 
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FIG. 10: Time dependence of the occupation of the exited level of the decaying system. Solid line shows results of numerical 
calculation, dashed line shows exponential decay according to Fermi's Golden Rule. The dotted line shows approximation 
according to Eq. H90H . The parameters used are the same as in Fig. [0] 
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600 800 
t/At 

FIG. 11: Typical quantum trajectory of the measured decaying system (solid line) when the detector interacts with the excited 
level. Figure shows the probability p ee that the measured atom is in the excited level. The dashed line shows exponential decay 
according to the Fermi's Golden Rule in the measurement-free evolution. The parameters used are the same as in Fig. 



Equation (|58l) in this case does not give correct decay rate of the measured system. In order to estimate the decay rate, 
we solve the Liouville-von Neumann equation ihp = [H, p] for the density matrix of the system with the Hamiltonian 
(|65f) - l|68|) . including additional terms describing decay of the non-diagonal elements with rate 1/tm, i.e., 



Pe0,e0 = -i^2(g(k)p gk ,e0 ~ Pe0,gkg(k)*), 
k 

Pgk,gk> = -i^kk'Pgk.gk' - i(9(k)*PeO,gk' - Pgk,e09(k')), 
PeO.gk = \-i(u)A-Wk) ) PeO.gk ~ i(5^ g{k')p g k',gk - Pe0,e0g(k)), 

Pgk,e0 = \-i{uk-U A ) ) Pgk.eO - i(g{k)*p e 0.e0 Pgk,gk>g(k')*). 



(94) 

(95) 
(96) 

(97) 



We solve equations (I94H - I|97|) using the Laplace transform method. Eliminating p e o.gk and p g k,e0 from the Laplace 
transform of Eqs. (|94|l - (|97|l one gets the equations for the Laplace transforms of the matrix elements of the density 
matrix, 



zpeo 



e0 (z)-l = -$>(A0| 2 ( 

k \ 



(z + iU) k k')pgk,gk'{z) 



1 



1 



z + i(u k - u>a) + 777 z + i(u A - uj k ) + ^ 
1 1 \ 



PeO,eo(z) 



k.k' 



z + i{LO k -uj A ) + z + i(w A - 0Jk') + 77 



TM 



xg(k)g(k')*p~gk,gk>{z), 
Y ( 9{k')g{k"Y 



— Pgk,gk"(z) + 



g(k")g(ky 



+g(k')g(ky 



1 



z + i(uj A - w k >) + — 
1 



— Pgk",gk'{z) 



Z + i(u k - LO A ) + — Z + i(LU A - U) k ') + — 



PeO,eo{z)- 



(98) 



(99) 
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FIG. 12: Time dependence of the occupation of the exited level of the decaying system when the interaction with the reservoir 
modes is described by Eq. 1781 . Solid line shows results of the numerical calculation, dashed line shows exponential decay of the 
measurement- free system with the decay rate given by Eq. I|79|l . The dotted line shows approximation according to Eq. <f lOip . 
In the calculations we used a = 2, while other parameters are the same as in Fig. |5] 

On the r.h.s. of Eq. I|99|l we will neglect the small terms not containing p e o e o(z). Expressing p g k,gk'{z) via p e o,eo{ z ) 
from Eq. I|99[) . substituting into Eq. I|98|l and replacing the sum over k by an integral we obtain 



1 



PeO,eo(z) 



= z+ I cLujG(uj) 



1 



1 



z + i(w - ua) + ^ z + i(u A -u) + ^ 
1 



duj J duj 1 G{uj)G{uj' 
z + i{uj - uj a ) + + z + i{uj A - w>) + 



z + i{u> — lu') 
1 



(100) 



The value of z at which the r.h.s of Eq. I|100|) is equal to zero gives the decay rate. Using the expression l|93l) for G(ui) 
and keeping only the first-order terms of the expansion into series of the powers of A -1 we get the measurement- 
modified decay rate 



e^-le^ 9 l 1 nA (5a 1) +l e ^^ ^ ■ 



(101) 



Equation l|101|) is valid only for sufficiently large duration of the measurement tm, since expansion into series requires 
that Atm ^S> 1. From Eq. I|101|) one can see that in order to obtain the quantum anti-Zeno effect the parameter a 
should be greater than 1. When the parameter a is less than 1 we get the Zcno effect, and when a — 1 the decay rate 
coincides with the decay rate of the free system. 

The probability that the measured system stays in the excited state, obtained from numerical simulation is presented 
in Fig. 1121 Figure clearly demonstrates the quantum anti-Zeno effect, the decay rate of the measured system is bigger 
than that of the measurement-free system. 
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IX. CONCLUSIONS 

We model the evolution of the measured quantum system as a detector using a two level system interacting with 
the environment. The influence of the environment is taken into account using quantum trajectory method. The 
quantum trajectories produced by stochastic simulations show the probabilistic behavior exhibiting the collapse of 
the wave-packet in the measured system, although the quantum jumps were performed only in the detector. Both 
quantum Zcno and anti-Zeno effects were demonstrated for the measured two level system and for the decaying system. 

The results of the numerical calculations are compared with the analytical expressions for the decay rate of the 
measured system. It is found that the general expression |J5SJ, obtained in Ref. [53], gives good agreement with the 
numerical data for the measured two level system and for the decaying one showing the quantum Zeno effect. Never- 
theless, when the interaction of the measured system with the reservoir is strongly mode-dependent, this expression 
does not give the correct decay rate. The decay rate in this case was estimated including additional terms describing 
decay of non-diagonal elements into the equation for the density matrix of the measured system and a good agreement 
with the numerical calculations is found. A good agreement of the numerical results with the analytical estimates of 
the decay rates of the measured system shows that the particular model of the detector is not important, since the 
decay rates mostly depend only on one parameter, i.e., the duration of the measurement l|22|). 
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